Transmissivity directional hysteresis of a nonlinear metamaterial slab with very small 

linear permittivity 
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We investigate propagation of a transverse magnetic field through a nonlinear metamaterial slab 
of sub-wavelength thickness and with a very small and negative linear dielectric permittivity. We 
prove that, for a given input intensity, the output intensity is a multi- valued function of the field 
incidence angle so that the transmissivity exhibits angular multi-stability and a pronounced direc- 
tional hysteresis behavior. The predicted directional hysteresis is a consequence of the fact that 
the linear and nonlinear contributions to the overall dielectric response can be comparable so that 
the electromagnetic matching conditions at the output slab boundary allow more than one field 
configurations within the slab to be compatible with the transmitted field. 



Optical hysteresis behavior and bistability are fasci- 
nating nonlinear phenomena which have attracted a large 
research interest in the last decades [l[ mainly for their 
potential photonic applications as optical memories, logic 
gates and optical computing devices [2j. The possibility 
of engineering the dielectric permittivity and the mag- 
netic permeability offered by metamaterials ^ has re- 
cently allowed to prove that the feedback mechanism sup- 
porting bistability is facilitated by the the opposite direc- 
tionality of the phase velocity and the energy flow in the 
negative index metamaterial Q. Nonlinear propagation 
of light through alternating slabs of positive and negative 
refractive index materials characterized by a vanishing 
average refractive index has been considered and bistable 
switching [5] and gap soliton formation [6] have been pre- 
dicted. Metal-dielectric multilayer structures have also 
proved to be media exhibiting interesting bistable be- 
havior [7] since the small average dielectric permittivity 
combined with the nonlinearity allows the sign of the ef- 
fective overall dielectric constant to be dependent on the 
optical intensity [H, ^ . 

In this Letter, we investigate electromagnetic trans- 
mission through a metamaterial slab of sub-wavelength 
thickness, characterized by a very small and negative lin- 
ear dielectric permittivity and exhibiting focusing Kerr 
nonlinearity. We numerically solve Maxwell equations for 
the problem of reflection and transmission of an inclined 
incident transverse magnetic (TM) plane wave and we 
show that the slab transmissivity is, for a given input 
intensity, a multi-valued function of the field incidence 
angle. The novel directional hysteresis behavior is due 
to the fact that the slab can host the extreme nonlin- 
ear regime where the linear and nonlinear contributions 
to the overall dielectric response can be comparable [ij 
so that the electromagnetic matching conditions at the 
output slab boundary allows more than one possible field 
configurations compatible with the transmitted field. 

Consider a transverse magnetic (TM) monochromatic 
field (with time dependence exp(— iuot)) E = E x (x, z)e x + 



FIG. 1: (Color online) (a) Geometry of the metamaterial slab, 
(b) Geometry of the TM incident (i), reflected (r) and trans- 
mitted (t) plane waves together with the TM field within the 
slab. 



E z (x, z)e z and H = H y (x,z)e y propagating through 
the metamaterial slab, embedded in vacuum, reported 
in Fig. 1(a) of sub- wavelength thickness (along the z- 
axis) L = 1/ko (fco = 00 /c) and consisting of alternat- 
ing, along the y-axis, two nonlinear isotropic and non- 
magnetic Kerr layers of thicknesses d\ and c?2- If the 
spatial period is much smaller than the vacuum wave- 
length (i.e. d\ + d2 <C 27r/fco) ? the slab behaves as an 
homogeneous medium characterized by the dielectric re- 
sponse (constitutive relations) 

{\E x \* + \E z \*)E x + \(El + El)El 

(\e x \ 2 + \e z \ 2 )e z + 1 -(eI + e*)e; 



D x = e ^eE x -\-x 
D z = eAeE ; 



X 



where the effective dielectric permittivity e = fe\ + (1 — 
f)e2 and effective nonlinear susceptibility \ — fXi + (1 — 
f)X2 are the averages of their layers counterparts (here 
/ = di/{d\ + cfo) is the volume filling fraction of layer 
1)@. Therefore, employing both nonlinear negative and 
positive dielectrics (i2e(ei) < and Refa) > 0), con- 
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ditions can be found so that the effective dielectric per- 
mittivity e has a negative real part (with a very small 
imaginary part) and very close to zero, and x > (see 
below for an explicit example) which is the situation we 
focus on in this Letter. Note that, due the the small 
value of e, the effective medium is able to support the 
extreme nonlinear regime where the nonlinear contribu- 
tions to Eqs.(pQ) can be comparable of even greater than 
the linear part eE [iQ. 

A TM plane wave (Ei,Hi) incoming from vacuum 
(z < 0) and impinging on the slab at z = with in- 
cidence angle 0, produces a TM reflected plane wave 
(E r ,H r ) for z < and a TM transmitted plane wave 
(Et,Ht) for z > L, as reported in Fig. 1(b). The fields of 
the three waves are E s = E s (e y x k s /fco) exp (zk s • r), 
H s = E s y/eo/voey exp (ik s • r) where s = (i,r,t), 
Ei, E r and E t are the amplitudes of the three 
waves and = = ko (sm0e x + cos0e z ), k r = 
ko(sin0e x — cos0e z ) are the wave vectors of the three 
waves. The TM electromagnetic field within the slab 
is of the form E = [E x (z)e x + E z (z)e z ] ex.p(ikox sin 0), 
H = [H y (z)e y ] ex.p(ikoxs'm6) so that, Maxwell equations 
V x E = zcj/ioH, V x H = —iuj'D yields the system 



dE x 

dz 
dH y 

dz 

(fco sm 0)H y 



= i(ko sm0)E z + iufioH y 
iuD x , 



(2) 



where D x and D z are given by Eqs.(pQ) and where the 
third equation implicitly defines E z as a function of E x 
and H y . In order to evaluate the slab transmissivity 
T = \E t \ 2 /\Ei\ 2 , we specify the amplitude of the trans- 
mitted field E t at the output face z = L and successively 
solve Eqs.(j2j) all the way to the input face at z = 0. At 
the slab input and output faces we require the continu- 
ity of the tangential components of the electric and mag- 
netic fields, namely (Ei — E r ) cos 6 = E x (0), yj eo//io(^i + 
E r ) = H y (0) (at z = 0) and E X (L) = E t cos#exp(z cos#), 
H y (L) = ^/eof/jLoEt exp(i cos#) (at z = L) where use has 
been made of the fact that L = 1/ko. In order to eval- 
uate the transmissivity for a given input intensity as a 
function of we have first evaluated a number of ampli- 
tudes Ei for various E t and and subsequently collected 
all the configurations corresponding to the same Ei . The 
resulting slab transmissivity T as a function of 6 is re- 
ported in Fig.2 for \Ei\ = 10~ 6 |x| — 1/2 (dashed line) and 
for |^| = 5.3 • 10- 4 |xr 1/2 ( solid line ), for a slab with 
e = —0.001 and \ > 0. The dashed line transmissivity 
is evaluated for an input intensity so low that it coin- 
cides with the transmissivity of the linear slab and it has 
been reported for comparison purposes. Note that the 
solid line transmissivity is a multi- valued function of the 
incident angle since it is threefold in the angular range 
-0.028 < 6 < 0.028 (radians). 
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FIG. 2: 

of 6 for \Ei\ = 10" 6 |xr 1/2 (dashed line) and for \Ei\ = 5.3 
10" 4 |xr 1/2 ( solid line ) 5 for a slab with 6 = -0.001 and x > 0. 



In order to physically grasp the predicted novel kind 
of multi-stability it is convenient to consider the field 
matching conditions at the slab boundaries. Exploit- 
ing the third of Eqs.((2]) it is evident that the continu- 
ity of H y across the slab boundaries is equivalent to 
the continuity of D Zl so that after setting E Z (L) = 
A z (E t /\E t \) exp(z cos#) (where the amplitude A z is real) 
and exploiting the second of Eqs.(pQ) together with the 
continuity of E x , the matching condition D Z (L~) — 
D Z (L+) yields 



6+^x(\E t \ 2 cos 2 + A 2 z 



A z + \Et\sm0 = 0, (3) 



which is a relation that, for a given E t and yields the 
possible values of A z . Note that, in the linear regime 
(X = 0) Eq.© yields A z = -\E t \ sinO/e i.e. the bound- 
ary conditions allow only one possible value of A z . The 
possible values of Ixl 1 ^ 2 ^ obtained by Eq.(j3]), in the 
case e = —0.001 and x > 0? are plotted in Fig. 3 from 
which we note that the surface is folded so that there 
are regions where, for a given pair Ixl 1 ^ 2 !^! and there 
are more than one accessible value of \\^l 2 A z . In other 
words, there are more than one possible electromagnetic 
fields compatible with the same output field and inci- 
dence angle. The origin of the surface folding is easily 
discussed by considering the situation of vanishing in- 
cidence angle (0 = 0) for which Eq.(|3|), in addition to 
the standard solution A^ = 0, also admits the solutions 
\X\ 1/2 A Z ±] = ±^-sign( X ) (2e/3 + xl^| 2 ) which are re- 
ported with two dotted lines in Fig. 3 on the upper and 
lower surface sheets. The surface is folded if Ixl 1 / 2 ^!^ 
is real (existence of the upper and lower surface sheets) 
and therefore a necessary condition is sign(ex) = — 1. 
It is worth noting that |x| 1/2 4 ±) are obtained by the 
vanishing of the term in square brackets of Eqs.(j3]) (for 
= 0) and this corresponds to requiring the nonlinear 
contribution to exactly compensate the linear part in the 
second of Eqs.(pQ). Therefore, as opposite to the situa- 
tion we are considering in this Letter where |e| <C 1, in 
standard materials where e is of the order of unity, the re- 
quired nonlinearity would be so large to effectively forbid 
the folding of the surface (and the consequent feasibility 



3 




FIG. 3: (Color online) Surface \x\ 1/2 A z obtained by Eq.©, 
in the case e = —0.001 and x > 0- 

of the directional multi-stability). Note now that each of 
the point of Fig. 2 corresponds to a field configuration sat- 
isfying the matching conditions so that the curve of Fig. 2 
has an image on the surface of Fig. 3 which we report as a 
solid line. The above mentioned multi-valuedness of the 
slab transmissivity of Fig. 2 is therefore explained since 
the different values of transmissivity at a give corre- 
sponds to different fields with the same input intensity 
and having different output longitudinal component A z 
belonging to different sheets of the folded surface. 

The multi-valuedness of the slab transmissivity is ac- 
companied by a directional hysteresis behavior. Suppose 
the incident plane wave to be initially normally imping- 
ing on the slab (0 = 0) so that the state A (with A z = 0) 
of Fig. (3) and Fig. (2) is excited since before the illumi- 
nation the slab was not polarized. By increasing the in- 
cidence angle 0, the electromagnetic state continuously 
varies toward the state B (see Fig. 3) and the transmissiv- 
ity becomes to depart from its linear behavior (see Fig. 2) 
since A z arrives at the upper sheet of the surface where 
its value is greater than its linear counterpart. From the 
state B, if is further increased, from Fig. 3 it is evident 
that the state has to undergo a sudden jump to the state 
C (belonging to the lower sheet of the surface) since there 
is no allowed state on the upper sheet which is continu- 
ously joined with B at a greater 0. Correspondingly, the 
transmissivity undergoes a sudden jump to a higher value 
(T(C) > T(B), see Fig. 2) since A z has a large magnitude 
on the lower surface sheet. Suppose now that, from the 
state C, is decreased so that the electromagnetic state 
now continuously varies toward the state D on the lower 
sheet (see Fig. 3). Correspondingly, the transmissivity 
varies along the curve from C to D of Fig. 2 assuming 
values different from those of the forward path from A to 
B, i.e. directional hysteresis behavior occurs. Evidently, 
the procedure can be extended to the states E-F-G-B and 



so on. Note that the three states A, D and G, while hav- 
ing different values of A z (see Fig. 3), share a common 
value of the transmissivity (see Fig. 2) since they have 
= so that, from the first of Eqs.(|2|), E z can not affect 
electromagnetic propagation. 

As a specific example, consider a TM field of wave- 
length A = 810 nm so that the slab thickness is L = 
1/fco = A/27T = 129 nm. We choose the slab layers of 
type 1 and 2 (see Fig. 1(a)) to be filled with silver (1) and 
AlGaAs (2), respectively, so that e\ = —28.81 + 1.55i, 
Xi = 3.46 x 10- 16 m 2 /V 2 and e 2 = 10.9 - 0.5865i, 
X2 = 6.56 x 10 -18 m 2 /V 2 . Note that we are exploit- 
ing the very large silver nonlinear susceptibility [7] and 
the fact that AlGaAs optically amplifies the radiation 
at the chosen wavelength if the sample is pumped by 
ultra-violet light so that the above negative value of 
Imfa) can be attained by adjusting the pump laser in- 
tensity [l0f| . Choosing d\ = 14 nm and d2 = 37 nm, we 
obtain e = -0.001 - O.OOOOlz and x = 10~ 16 m 2 /V 2 
which are values adequate to observe the directional 
hysteresis. The intensity of the incident plane wave is 
/ = (1/2) ^ 6q / no\Ei\ 2 which, for the amplitude E{ cho- 
sen for evaluating the transmissivity of Fig. 2 and the 
above effective nonlinear susceptibility x, yields the value 
/ = 372 W/cm 2 , an extremely small value if compared to 
those (~ MW/cm 2 ) normally required for observing the 
standard optical bistability. 

In conclusion, we have shown that a sub-wavelength 
nonlinear metamaterial slab with a very small linear per- 
mittivity is characterized by a transmissivity exhibiting 
a novel kind of directional hysteresis behavior which is 
observable at very low intensities of the incident field. 
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